Multi-band superconductivity in LaFeAsOo.gFo.i single crystals probed by high-field 

vortex torque magnetometry 
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To probe manifestations of multiband superconductivity in oxypnictides, we measured the angular 
dependence of the magnetic torque t(6) in the mixed state of LaOo.gFo.iFeAs single crystals as a 
function of temperature T and magnetic fields H up to 18 T. The paramagnetic contribution of the Fe 



ions is properly treated in order to extract the effective mass anisotropy parameter 7 = (m c /m a i 



u ) 1/2 

from t{9). We show that 7 depends strongly on both T and H, reaching a maximum value of ~ 10 
followed by a decrease towards values close to 1 as T is lowered. The observed field dependencies of 
the London penetration depth \ a b and 7 suggest the onset of suppression of a superconducing gap 
at H « H c2 /3. 



PACS numbers: 74.25.-q, 74.25. Ha, 74.25.Op, 74.70.Dd 
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I. INTRODUCTION 

The recently discovered superconducting 
oxypnictides^ have similarities with the high T c 
cuprates, such as the emergence of superconductivity 
upon doping a parent antiferromagnetic compound.— ~— 
Theoretical models^ suggest unconventional pairing, 
consistent with Andrcev spectroscopy^ penetration 
depths and photocmission measurements^ which in- 
dicate nodless s-wave pairing symmetry. Experiments 
at high magnetic fields^r— penetration depth£ ' n i 13 
and heat capacity measurements^ are consistent with 
multiband scenarion^ Other models suggest an 
important role for electronic correlations^ and even a 
possibility of unconventional pairing mechanisms*^ 

The comparatively high critical temperatures T c 
and extremely high upper critical fields H C 2 of the 
oxypnictides^^ indicate their promising prospects for 
applications if, unlike the layered cuprates, a sizeable vor- 
tex liquid region does not dominate their temperature- 
magnetic field (T — H) phase diagram. It is therefore 
important to reveal the true behavior of the anisotropic 
magnetization in the vortex state of the oxypnictides, 
particularly the extent to which vortex properties are af- 
fected by strong magnetic correlations and multiband ef- 
fects. For instance, multiband effects in MgB2 can mani- 
fest themselves in strong temperature and field depen- 
dencies for the mass anisotropy parameter 7(T, H) = 
{riic/mai,) 1 / 2 and the London penetration depth \{T,H) 
even at H <C -ff C 2 ' 17lls Similar effects in pnictides would 
be consistent with the multiband pairing scenarios^ Yet, 
there are significant differences between two-band super- 
conductivity in MgB2 and in oxypnictides: in MgB2 the 
interband coupling is weak, while in the oxypnictides it is 
the strong interband coupling which is expected to result 
in the high T c Jl Thus, probing multiband superconduc- 
tivity in oxypnictides by magnetization measurements 
requires high magnetic fields, which can suppress the 
supcrfluid densities in both bands by circulating vortex 



currents. In this work we address these issues, present- 
ing high-field torque measurements of the anisotropic re- 
versible magnetization in LaFeAsOo.gFo.i single crystals. 
Our measurements of j(T, H) up to 18T and extended 
temperature range, 4 < T < 15 K reveals a different be- 
havior in j(T,H) as c omp ared to recent low- field torque 
measurements^ Ref. [UJ shows a 7 that increases con- 
tinuously as T is lowered reaching a maximum value of 
~ 20 for both NdFeAsO . 8 F .2 and SmFeAs0 . 8 Fo.2 sin- 
gle crystals. In contrast, our results indicate that 7 in 
LaFeAsOo.gFo.i reaches a maximum of ~ 10 decreasing 
asymptotically towards 1 as T is lowered. 

Measurements of the equilibrium magnetization 
m(T, H) of the vortex lattice in LaOo.gFo.iFeAs are com- 
plicated by the smallness of m(T, H) caused by the large 
Ginzburg-Landau parameter, k = A/£ > 100 and by 
the background paramagnetism of the normal stated 
which can mask the true behavior of m(T,H). In this 
case torque magnetometry is the most sensitive technique 
to measure the fundamental anisotropy of the param- 
eters of nt{T,~3) in small single crystals. The torque 

"7* = Horrl X ~H acting upon a uniaxial superconductor is 
given by 



t{6) 



HV(f> {7 2 ~ 1) sin 26> 
167r M oA 2 b 7 1/3 £(#) 
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where V is the sample volume, 4>q is the flux quan- 
tum, is the upper critical field along the ab planes, 
r\ ~ 1 accounts for the structure of the vortex core, 
is the angle between if and the c-axis, e{9) = (sin 2 9 + 
7 2 cos 2 6 1 ) 1 / 2 and 7 = A C /A Q (, is the ratio of the London 
penetration depths along the c-axis and the ab-plane. 
The first term in Eq. ([!} was derived by Kogan in the 
London approximation valid at H c i <C H <?C -ff C 2,— while 
the last term describes the torque due to the background 
paramagnetism, for which r m = fio(Xc — Xa)VH 2 /2 and 
Xc and Xa are the normal state magnetic susceptibilities 
along the c-axis and ab- plane, respectively In general, 
7 = (mc/iriab) 1 K/Kb, (where m c and m ab are the 



2 



Ginzburg-Landau superconducting effective masses) but 
both ratios are assumed to be equal in the model lead- 
ing to Eq. [TJ As will be shown below, the paramagnetic 
term in Eq. ^ in LaOo.gFo.iFeAs can be larger than 
the superconducting torque, which makes extraction of 
the equilibrium vortex magnetization nontrivial. In this 
work we propose a method, which enables us to resolve 
this problem and measure the true angular dependence 
of the superconducting torque as a function of both H 
and T, probing the concomitant behavior of 7(T, H) and 
X a b(T,H) and manifestations of multiband effects. 



II. EXPERIMENTAL 

Undcrdoped single crystals of LaOi-zF^FcAs with 
typical sizes of 80 x 60 x 5 urn 3 were grown by the flux 
method described in Ref. [2l|. The samples had a criti- 
cal temperature T c ~ 15 K as determined by the SQUID 
magnetometry, and as shown in Fig. 1. The width of 
superconducting transition, measured by a commercial 
SQUID magnetometer under a field H = 10 Oe after 
cooling the crystal under zero field, is AT C ~ 3.5 K. This 
relatively broad transition may not reflect the sample 
quality but mostly results from the penetration of vor- 
tices in a plate-like crystal which has a large demagne- 
tization factor and thus reduced lower critical field. Al- 
though in some crystals from the same batch the width 
of the resistive transition, from the very onset of the re- 
sistive transition to the zero resistance state, is observed 
to be as large as AT ~ 3 to 4 K , consistent with the 
values reported in the literature for crystals with simi- 
lar composition^. The fraction of F quoted here, corre- 
sponds to a nominal value since its precise content is very 
difficult to determine in such small single crystals. How- 
ever, a superconducting transition temperature T c ~ 15 
K, see Fig. 1, firmly places these crystals within the un- 
dcrdoped state, following the overall phase diagram dis- 
playing T c (x) as a function of the F content As ar- 
gued in Ref. F is expected to be inhomogeneously dis- 
tributed throughout the samples. Samples were attached 
to the tip of a piezo-resistive micro-cantilever placed in a 
rotator inserted into a 3 He cryostat. The ensemble was 
placed into a 18 T superconducting solenoid. Changes 
in the resistance of the micro-cantilever associated with 
its deflection and thus a finite magnetic torque r was 
measured via a Wheatstone resistance bridge. 



III. RESULTS AND DISCUSSION 

Fig. 2 (a) shows the angular dependence of the torque 
t{9) at T = 8 K and under H = 5 T for a LaFeAsO . 9 F .i 
single crystal. Hysteresis, resulting from the irreversible 
magnetization due to vortex pinning is observed be- 
tween increasing and decreasing angular sweeps. Black 
line depicts the average value of both traces, T rev (6) = 
(t(0)up + T"(^)down)/2 defined as an equilibrium magnc- 
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FIG. 1. (color online) (a) Magnetization M as a function of 
temperature for a typical LaOo.gFo.iFeAs single crystal mea- 
sured under an external field of 10 Oe. The small mass of the 
single crystal, of just a few micrograms, leads to the observed 
large scattering in the data points. Red markers correspond 
to M(T) averaged over 20 raw data points. 



tization, where 9 up and #down indicate either increasing 
or decreasing angle sweeps, respectively. This reversible 
torque contains both superconducting and paramagnetic 
contributions. Given that vortex pinning and concomi- 
tant hysteresis are bound to disturb the equilibrium mag- 
netization, we have checked that our definition of t icv (9) 
leads to reproducible traces, by re-measuring t{9) for 
both 9 up and #down within two distinct angular ranges, 
as shown in Figs. 2 (b) and 2 (c). The respective T lev (9) 
traces, i.e. blue and magenta traces, are plotted in Fig. 
2 (d) together with the original T rcv ((9) from Fig. 2 (a). 
As seen, all traces overlap almost perfectly. Both sharp 
spikes seen at either side of = 90° are likely due to 
pinning effects. 

Fig. 3 (a) shows the angular dependence of the 
torque t{9) at T = 10 K and H = 5 T for a sec- 
ond LaFeAsOo.gFo.i single crystal. In this crystal and 
at this temperature one does not observe any structure 
that might be attributable to pinning and which would 
compromise a fit attempt of the resulting r rov ((9) (brown 
line) to either Eq. (1) or Eq. (2). As shown in Fig. 3(b), 
the paramagnetic term r m sin 29 increases rapidly as the 
field increases, preventing a direct fit of T rQV (0) to Eq. 
(1) since r m sin2# interferes with the sin 29 harmonics 
of the first term. Yet the superconducting parameters 
7(T, H) and X(T,H) can be unambiguously extracted 
from the data adding a 90° translation of r lev (9) to itself, 
i.e. T rO v(0) + t tov (9 + 90°), where the paramagnetic term 
in t(9) cancels out: 



t{9) + t(9 + 90°) 



U^o(7 2 - 1)# sin 20 



1 ln (vHi\ * 



e(0) \s(9)H J e*(0) 



le^oA^ 1 / 3 
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FIG. 2. (color online) (a) Magnetic torque r for a LaOo.gFo.iFeAs single crystal for increasing and decreasing angle (9) sweeps 
(red lines), at H = 5 T and T = 8 K. Arrows indicate either increasing (9 up ) or decreasing (Odown) angles. Black line is the 
reversible torque component or T rB v(6) defined here as the average between both traces, (b) Same as in (a) but in a smaller 
angular range, (c) Same as in (a) and in (b) but in an even shorter angular range, (d) The resulting r lev (9) from Figs, (a), (b) 
and (c). 



where e*(6) = (cos 2 9 + 7 2 sin 2 6*) 1 / 2 . This procedure is 
illustrated by Fig. 3 (c) which shows r + (9) as a function 
of 9 for all curves in Fig. 3 (b) . Notice that the amplitude 
of t + (9) is considerable smaller than that of r lev (9). 

Since t+(0) in Eq. (2) depends on three fit param- 
eters, 7, A and rjH C 2, different sets of parameters may 
give equally good descriptions of the experimental data. 
One can circumvent this difficulty by extracting r]H C 2 
from the amplitude of t + (9)/H plotted as a function 
of ln(l/7J) for several values of 9, as shown in Fig. 
4 (a). If at a given temperature, these measurements 
are performed up to high enough fields, the extrapo- 
lation of t + (9)/H to zero evaluates the value H+ at 
which (nH^ /e(6)Hj\ = 1. From the two extrapolated 
values of H* for 9 X = 110° and 9 2 = 120°, we ex- 
clude 7 and obtain (r/iJ c2 ) 2 = H 2 (9 1 )H 2 (9 2 ) sin(6> 2 - 
0i)sin(0 2 + 9 l )/[Hl{9 l )cos 2 9 1 - H 2 (6 2 ) cos 2 6>i], which 
yields rjH^ ~ 27 T for T = 10 K. This simple method 
provides a thermodynamic estimate for H^ 2 without the 
need of extremely high magnetic fields, H ~ H C 2(9). Fur- 
thermore, with r]H^2 — 27 T we obtain excellent and sta- 
ble two-parameter fits of t+(6) to Eq. (2) for virtually all 
field values, as shown in Fig. 4 (b). Here the parameters 
7 and A no longer interfere in the fit as A affects only 
the magnitude but not the shape of r + (9) . We confirmed 



that the extracted parameters lead to a small difference 
in the logarithmic terms of Eq. (2) at H = H+. 

Using the method outlined above, we extracted the 
field and the temperature dependencies of j(T 7 H), 
\ a b{T,H). An example, is given by the orange line in 
Fig. 5 (a) for H = 6 T and T = 10 K. Having fixed 
7, A, and r/H^ we can now fit the original t icv (9) to 
Eq. (1) (red line), leaving the amplitude of the param- 
agnetic component t,„ as the only adjustable parameter 
(magenta line). Figure 5 (b) depicts the resulting ampli- 
tudes for both r rn and the superconducting contribution 
tk (x H/ A 2 b in Eq. (1) as a function of field for T = 10 K. 
The amplitude of tk follows the expected linear in field 
dependence up to H = 8 T, from which point it starts to 
decline continuously. On the other hand, r m follows the 
expected H 2 dependence characteristic of the torque of 
an anisotropic paramagnetic background. To our knowl- 
edge, this is the first magnetometry method which allows 
an independent extraction of the parameters of a mag- 
netic superconductor. 

The curves t{9) measured for T > 4 K, indicate 
that riH^{T) is described by r)H$(T) = r/iJ c a 2 b (0)[l - 
(T/T c ) 2 ]. As seen in Fig. 5, the irreversible component 
in t(9) grows quickly as the temperature is lowered. At 
T = 4 K despite the observed large irreversibility, the 
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FIG. 3. (color online) (a) Magnetic torque r for a 
LaOo.gFo.iFeAs single crystal for increasing (clear blue line) 
and decreasing (blue line) angle (0) sweeps, at H — 3 T and 
T = 27 K. (b) T rcv (6») for several field values at T = 10 K. 
(c) The symmetrized component r lev (6) + r lev (8 + 90°) of the 
torque due solely to the reversible vortex magnetization. 
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FIG. 4. (color online) (a) The amplitude of [r lev (8) + T Tev (9 + 
90°)]/H as a function of ln(H~ 1 ) for several angles, taken 
from the data in Fig. 1 (c). Red lines are linear extrapolations 
of [Trev(0) + r re v(0 + 90°)]/ H -»■ for 9 = 110° and 120°, 



which yields rjH" 



27 T for LaFeAsOo.gFo.i at T 



10 



K. (b) With the value ■qH'g ~ 27 we obtain excellent fits of 
r rev (0) + t icv (8 + 90°) to Eq. (2) for all field values. 



resulting reversible component in t{6) is still nearly per- 
fectly described by Eq. (1), see Fig. 6 (a). However, as 
the temperature is lowered to T = 1.5 K additional struc- 
tures emerge in t(9) for 9 close to 90°, see Fig. 6 (b). We 
ascribe these features to the intrinsic pinning of vortices 
by the planar structure of the material, which prevents 
a reliable extraction of the reversible component in t{9) 
below T ~ 4 K. A detailed analysis and discussion will 
be provided elsewhere^ 

Figure 7 (a) shows the resulting temperature depen- 
dence of the total superfluid density n s oc A~ b 2 mea- 
sured under a field H = 5 T. A field of 5 T is high 
enough to suppress the irreversible component in the 
torque as compared to the reversible one, and to ne- 
glect demagnetization factors and geometrical barriers 
for the penetration of magnetic flux through the sam- 
ple. Our torque magnetometer was not calibrated, so 
we could only measure the temperature dependency of 
^ab{T) but not the absolute values of X a b or n s . Despite 
the unavailability of points at the lowest Ts, the tem- 
perature dependence of n s (T) exhibits a positive curva- 
ture as seen in penetration depth measurements in the 
compounds i?FcAsOo.gF .i (iZ=La, Nd, with T c s of 14 
and 45 K, respectively)^ But it contrasts with penetra- 



tion depth results for PrFeAsOi-^ (T c ~ 35 K)2£ and 
for SmFcAs0 .8Fo.2 (T c ~ 44) YX, which finds evidence 
for a superconducting state characterized by more than 
one (non-nodal) gap. In Ref. [26| such an upward cur- 
vature was explained in terms of the s ± scenario, with 
the superconducting gap ratio A1/A2 — 1/3. As fol- 
lows from Fig. 7 (b), the associated A a &(T) in our lim- 
ited temperature range follows a T 2 dependence at low 
and intermediate temperatures, in contrast to the lin- 
ear dependence expected for a nodal gap as observed 
in LaFeP02£, or the nearly exponential behavior ex- 
pected for a clean s-wave superconductor, as reported for 
SmFeAsOo.sFo.^ and PrFeAsOi_ y i2I One important de- 
bate concerning the Fe arsenides is precisely the origin of 
these differences, since a common pairing mechanism^^ 
has been proposed for all compounds based on the pre- 
sumed similarity of their electronic structure. At the 
same time, it was recently argued that the extended s- 
wave scenario can lead to either nodeless or nodal gaps, 
depending on the interplay between intraband and in- 
terband interactions.— Penetration depth measurements 
over an extended region in reduced temperatures (T/T c ) 
also revealed a Independence for A a b in i?FcAsOo.gF .i 
(R =La, Nd)2£ as well as in Bai_ 2; K 2; Fe2As2 3 This 
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FIG. 5. (color online) (a) Reversible torque r lcv (8) for a 
LaFeAsOo.gFo.i single crystal (blue markers) at H — 6 T 
and T — 10 K. Orange line shows the reversible vortex torque 
tk and the magenta line shows the paramagnetic torque ex- 
tracted from the data shown in Fig. 3 using the deconvolution 
procedure described in the text. Red line shows the sum of 
tk with the paramagnetic term, (b) The amplitude of tk 
normalized by (-y 2 — 1)/7 1//3 (blue markers) and that of the 
paramagnetic term (clear blue markers) as functions of H. 
At lower fields, tk(H,8) increases linearly with H, while the 
paramagnetic term exhibits a H 2 dependence. 



power-law temperature dependence was deemed to be 
consistent with the s± muti-gap scenario, if either strong 
interband impurity scattering^ or pair-breaking effects^ 
are important. Our observation of a n s oc T 2 is limited by 
the restricted temperature range imposed by the pinning 
effects at lower Ts, yet we believe that it is important 
to expose the agreement for the anomalous temperature 
dependence of the penetration depth between magnetic 
torque and surface impedance measurements in the case 
of LaFeAsOo.gFo.i. Given the width of the superconduct- 
ing transition ~ 3— 4 K, it is possible that local variations 
of T c due to an inhomogeneous distribution of F, could 
affect the dependence of n s on T at higher temperatures. 
In addition, given that our measurements of n s (T) were 
performed under a field of 5 T, we cannot completely rule 
out that the temperature dependence of H C 2(T) might 
influence the power law dependence for n s (T) reported 
here, particularly at higher temperatures. 

Figure 7 (b) also displays the temperature dependence 
of the mass anisotropy parameter 7, which starts at 7 ~ 7 
for T < T c , increases toward a maximum of ~ 10 at T ~ 
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FIG. 6. (color online) (a) r as a function of the angle 6 and 
respectively for increasing (blue line) and decreasing (orange 
line) angle scans at T = 4 K and under a field H = 5 T. 
Black line corresponds to the average between both traces, 
(b) Same as in (a) but for T = 1.5 K. Notice the emergence 
of sharp peaks in t(6) (red arrows) for 6 close to 90° resulting 
from the intrinsic pinning by the planar structure. 



0.75 T c , and then decreases to « 1 as T is lowered (as it 
happens in nearly all Fe based superconductors, see Refs. 
and [Lsh . Our unpublished transport measurements in 
LaFei-zCo^AsO indicate a very similar T-dependence 
for lH =H$/m 2 . 

Figures 7 (c) and 7 (d) show the field dependencies of 
7 and n s for T = 10 K. Here 7 increases by nearly a 
factor of 2, from ~ 8 at small fields to^l5atif = 18T. 
On the other hand, n s remains nearly constant in fields 
up to H v = 8 T, and then decreases at higher fields, ex- 
trapolating to n s = at .ff (10 K) w 25 T. The complete 
suppression of n s at this value is consistent with our pre- 
vious estimate of r]H C 2 « 27 T obtained above. The sup- 
pression of n s (H) upon application of magnetic field is 
consistent with the onset of orbital pair-breaking by cir- 
culating vortex currents as the spacing between the vor- 
tex cores becomes smaller than « y/3 of the vortex spac- 
ing at H C 2 where superconductivity is fully suppressed. 
For instance, the single-band Ginzburg-Landau theory 
predicts A~ 2 oc n s (H) ~ n s (0)(l — H/H c2 ) for T close to 
T c i22 Interestingly, the ratio H V /H C 2 ~ 1/3 turned out to 
be approximately equal to the ratio of the superconduct- 
ing gaps on the electron and hole pockets of the Fermi 
surface as observed in Ref. [26|. Around H = 6 - 7 T, 
there is a spike in A and/or a dip in n s , the origin of 
which remains unclear. 
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FIG. 7. (color online) (a) Temperature dependence of the 
superfluid density n s oc A~ 6 2 extracted from the torque data 
at H — 5 T. (b) Temperature dependence of both 7 (clear 
blue markers) and \ a t (blue markers). Here X a b(T) exhibits 
a quadratic dependence on T. (c) 7 as a function of magnetic 
field as extracted from the fittings in Fig. 2 (b). (d) The 
superfluid density n s oc X~? as a function of H at T = 10 K. 



We described our data using the simplest Eq. (1) for 
the torque in a uniaxial superconductor, assuming that 
the angular dependencies of A and £ are controlled by 
the single anisotropy parameter 7. A more complicated 
expression for r with two different anisotropy parame- 
ters -fn = H^/H^ 2 and 7^ = \ a b/X c was suggested 
for multiband superconductors^ However given that the 
formula for r of Ref. was obtained on phenomenolog- 
ical grounds, and the angular dependencies of A(#) and 
H C 2{9) derived from the multiband BCS theory are much 



more complicated j& we believe that the use of that ex- 
pression for t with different 7^ and 7^ may not unam- 
biguously reveal any new physics as compared to the sim- 
pler Eq. ([1]). Indeed, the fit of our data with r from Ref. 
|33| gave behaviors qualitatively similar to those shown 
above: 7# ~ 10 \t^t c increases to a maximum value 
of 18 before decreasing to much lower values as T — >• 0. 
Meanwhile 7a = 1.4 \tkt c exhibits a mild decrease to- 
wards 1 as T — > 0. This common tendency of 7 to ap- 
proach a value close to 1 upon decreasing T, observed on 
many Fe-based superconductors, indicates the increasing 
contribution of the Pauli pair-breaking effects at lower 
Ts which is neither taken into account in Eq. (1) nor in 
the expression for r of Ref. [33| . 

In conclusion, we report the torque measurements of 
reversible magnetization of LaOo.gFo.iFeAs single crys- 
tals from which the field dependent mass anisotropy pa- 
rameter 7 and the reduced London penetration depth 
\(H,T)/\(Q,T) were measured. The significant field de- 
pendencies of 7 and X(H,T)/X(0,T) are consistent with 
multiband pairing while the observed quadratic temper- 
ature dependence of X(H, T)/A(0, T) at low and interme- 
diate T is consistent with the pairing with interband 
impurity scattering. 
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